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According to the general classification of topological insulators, there exist one- dimensional chi- 
rally (sublattice) symmetric systems that can support any number of topological phases. We intro- 
duce a zigzag fermion chain with spin-orbit coupling in magnetic field and identify three distinct 
topological phases. Zero-mode excitations, localized at the phase boundaries, are fractionalized: 
two of the phase boundaries support ±e/2 charge states while one of the boundaries support ±e 
and neutral excitations. In addition, a finite chain exhibits ±e/2 edge states for two of the three 
phases. We explain how the studied system generalizes the Peierls-distorted polyacetylene model 
and discuss possible realizations in atomic chains and quantum spin Hall wires. 
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Introduction- The discovery of quantum spin Hall ef- 
fect in two dimensional electron systems initiated an 
enormous wave of research exploring topological prop- 
erties of gapped band structures [1, 2]. This has lead to 
various recent discoveries of novel states of matter, com- 
monly known as topological insulators, exhibiting strik- 
ing properties such as exotic electromagnetic responses 
and gapless surface states [3]. The quest for novel topo- 
logical insulators, guided by the comprehensive classifica- 
tion presented in Ref. [4] , is one of the most active topics 
in condensed matter physics at the moment. 

In one spatial dimension (ID) the recent attention has 
concentrated on the topological superconducting phases 
realized in spin-orbit coupled semiconductor nanowires 
in proximity of a supeconductor [5|, l6|] . It was proposed 
that networks of wires could serve as a platform for topo- 
logical quantum computation [7]. However, the study of 
topological phases in ID systems can be traced back over 
30 years to the discovery of soliton excitations in poly- 
acetylene [Sj. Indeed, the Su-Schrieffer-Heeger model is 
an example of an insulator exhibiting a chiral symmetry 
[9|. The zero-energy states located between differently 
dimerized domains support neutral spin- 1/2 excitations 
and spinless ±e excitations. 

In this work we introduce a zigzag fermion chain model 
with a spin-orbit coupling in a magnetic field. This sys- 
tem belongs to the BDI class in the classification scheme 
and has three distinct topological phases depending on 
the orientation and magnitude of the magnetic field. The 
domain wall (DW) excitations localized at the phase 
boundaries are fractionalized: two of the DWs support a 
single zero mode with possible charge states ±e/2 while 
one DW supports two zero modes leading to ±e and 
neutral excitations. Furthermore, a finite chain exhibits 
±e/2 charged edge states for two of the three phases. The 
most striking differences to the polyacetylene model in- 
clude existence of three distinct phases (instead of two) 
and possibility to realize ±e/2 excitations that are ab- 
sent in polyacetylene due to spin degeneracy. Zigzag 



chains are realized in a number of systems which indicates 
that the studied model could be relevant, for example, in 
monoatomic gold chains in magnetic fields. Additionally, 
a continuum limit of the considered model describes nar- 
row wires of quantum spin Hall insulators. Interestingly, 
in this case the ±e/2 excitations can be thought of as 
a bound pairs of ±e/4 charges on opposite edges of the 
wire. 

ID models and phase diagrams- In strictly ID chains 
the Peierls instability commonly results in geometrical 
deformations doubling the unit cell and creating a sublat- 
tice structure. If atoms are not restricted to a line, chains 
can relax to planar zigzag or more complicated struc- 
tures. Motivated by this, we consider a zigzag fermion 
chain shown in Fig. [Tal The Hamiltonian of the system 
is 

H = tY,clcj^S ^(-1)\^ + 
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where the first term describes the nearest-neighbor hop- 
ping and the second term encodes the hopping modu- 
lation due to dimerization. The third term arises from 
the next-nearest neighbor spin-orbit hopping and the last 
term corresponds to the Zeeman splitting m = (m^^, rriy) 
due to the external magnetic field. The spin-orbit term 
depends on the orientation of the two bonds connect- 



ing the next-nearest neighbors e^ 



di X dn and van- 



ishes for strictly linear chains [1]. The spin-orbit coupling 
and magnetization are crucial for the studied topological 
properties of the system while the hopping modulation is 
included for the sake of generality. The model ([1]), with 
an in-plane magnetic field leads to the Bloch Hamiltonian 



H{k) = 2tcos{ka)ax + 2S sm{ka)ay-\- 
-\-2Xsm{2ka)o^zSz + TnSx 
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Figure 1: (a) Studied ID models, a zigzag fermion chain with a spin-orbit hopping and narrow quantum spin Hah 

wire with a finite inter-edge couphng. Topological properties of these systems are closely related, (b-d) The phase 

diagrams of Hamiltonians ([2]) and (dj) as a function of model parameters, (b) For vanishing hopping modulation 

(5 = there is a (1/ — 1) interface at r?i = 0. (c) For a finite 5^ a trivial phase emerges for \m\ < \S\. (d) The 

Hamiltonian (jlj) shows three different phases for finite inter-edge coupling A. The normalization of (|3]) is multiplied 

by 2 since continuum models produce half-integer values. 



where Pauli matrices ai and Si operate in the sub- 
lattice and spin space, respectively. The Hamiltonian 
([2]) possesses a chiral symmetry: it anticommutes with 
C = (JzSy. This property has important consequences 
on the energy spectrum and topological properties of 
the system y. If H has an eigenstate \E^) with en- 
ergy Ei, there exists another eigenstate C\E^) belonging 
to the energy — £^i, so the spectrum is symmetric about 
E = 0. Zero modes, if existing, can be chosen to coin- 
cide with the eigenstates of the generator of the chiral 
symmetry, C\ip^) = :^\ip^). In odd spatial dimensions 
systems with a chiral symmetry may posses topologically 
nontrivial phases characterized by a Z- valued topological 
index introduced in Ref. 14]. Here we employ alternative 
but equivalent formulation [10|. In ID the invariant takes 
the form 



i^[H] = ^fdkTT [CH-^dkH] , 



(3) 



where the integration extends over the Brillouin zone for 
Bloch Hamiltonians and over the real line for continuum 
models. 

In Fig. [T]we plot the invariant ([3]) for Hamiltonian ([2]) 
as a function of the magnetic field, revealing the phase di- 
agram. The system has three distinct phases v = ±.1 and 
z/ = 0, separated by quantum phase transitions. In these 
interfaces where the invariant changes its value, there ex- 
ist zero energy states protected by the chiral symmetry. 
The phase of the system depends on m and can thus be 
changed by tuning the magnitude of the external field, 
i.e., by introducing magnetic DWs into the system. At 
large field values (|m| > 2t) the magnetization dominates 
and V = {). We will show below that this phase can be 
identified with vacuum. 

At vanishing spin-orbit coupling A = and magneti- 
zation m = 0, model (|2]) describes a Peierls-distorted 
polyacetylene chain [8|. The resulting Hamiltonian an- 
ticommutes with C = (Tz which serves as the generator 



of the chiral symmetry. Evaluation of the invariant ([3]) 
reveals two distinct topological phases with v = sign J. 
These phases correspond to the well-known degenerate 
dimer configurations of the chain. There is, however, im- 
portant differences between the polyacetylene model and 
the zigzag chain. In polyacetylene there is no z^ = 
phase which in our model is identified with vacuum. As 
a consequence, a finite polyacetylene chain does not ex- 
hibit robust edge states unlike v = ±.1 zigzag chains. 
In addition, due to the spin degeneracy in polyacetylene, 
there are always two zero modes at phase boundaries. As 
we emphasize below, due to this fact the zigzag chain ex- 
hibits fractional domain wall charges ±e/2 whereas poly- 
acetylene has ±e and neutral excitations. 

There is an interesting connection between the model 
([2j) and quantum spin Hall (QSH) wires (Fig. [Taj). A 
QSH bar supports helical edge states where the velocity 
and spin of the edge carriers are completely correlated 
[3] . In a narrow enough wire the edge modes localized on 
opposite edges begin to overlap, enabling interedge hy- 
bridization which gaps the edge spectrum. The Hamilto- 
nian describing the systems is 



H{k) = vpkazSz + Aax + mSx 



(4) 



where the first term describes the unperturbed edge 
states, the second term comes from the inter-edge tunnel- 
ing and the last term corresponds to magnetization. The 
Pauli matrices a now operate in the subspace of the two 
edges of the wire. Setting 6 = and expanding the band 
Hamiltonian ([2]) around A: = or /c = tt, and identifying 
2t ^ A and 4Aa -^ vp, one recovers (J4|). The close con- 
nection between (J2]) and (j4]) suggests that the topological 
properties of these systems bear similarities. The phase 
diagram of (|4|), shown in Fig. [id] reveals three distinct 
phases as a function of the relevant control parameters 
A and m. The zero-energy states at DWs are also very 
similar to those of the zigzag chain as we explain below. 
Zero modes and fractional excitations- Now we analyze 



zero modes by numerical diagonalization of the lattice 
model and analytical methods and identify associated 
quantum numbers. As mentioned above the spectrum 
is symmetric with respect to £^ = 0. This property and 
the number of the zero modes at DWs are sufficient to 
determine the charge states of the zero mode excitations 
[8|,lll]. According to the general argument, one-half of an 
electron state per zero mode is missing from the valence 
band in the vicinity of a DW. If a DW supports a single 
zero mode, then the possible charge states at half filling 
are ±e/2 depending on the population of the zero mode. 
If there exist two zero modes, possible charge states are 
±e and 0. The total energy of the system is independent 
of the population of the zero modes. In finite systems 
where the DWs come in pairs zero modes are split to 
linear combinations of \i^^) with finite energies around 
^ = 0, approaching zero as the separation of the DWs 
increases. 

At the DWs where the phase u = 1 changes to i^ = 
(referred as the 1/0 DW) we solve the model ([T]) in a finite 
ring with spatially varying magnetization m. For simplic- 
ity we assume that (5 = in this section, so there are two 
phase boundaries located at places where m crosses 2t. 
We choose a magnetization profile so that the DWs occur 
at the opposite sides of the ring. As illustrated in Fig. [2a] 
we discover two mid-gap energy states, both located at 
the DWs. When the ring size is increased the energies 
move closer to zero and eventually become degenerate. 
This result can be complemented by an analytical calcu- 
lation based on the continuum approximation of Eq. (|2]) . 
Linearizing ([2]) at Brillouin zone point k = (or tt) where 
the gap closes, we solve the Schrodinger equation for zero 
energy states. We find one normalizable solution V^i/o(^) 
('00/1 (^)) located at a DW (anti-DW), 
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where |a±) = (| t)|^±) + | ;)|^^)) and | t,;), \y±) denote 
the eigenstates of az and Sy. The DW profile is deter- 
mined by m{x). For '^i/q (V^o/i) '^{x) is an arbitrary 
decreasing (increasing) non- vanishing function satisfying 
m(xo) = — 2t and approaching constant value away from 
xq. The normalization constant A is determined by the 
detailed form of m{x). According to the general argu- 
ment, there is a charge aggregate of ±e/2 localized at 
the DW depending on whether the zero mode is popu- 
lated (— e/2) or not (e/2) [8]. Since the total number 
of electrons must be an integer, the DWs must come in 
"soliton-antisoliton" pairs which is clear in the ring ge- 
ometry. In a finite linear chain in the phase v = ±.1 ter- 
minated by hard-wall boundary conditions we discover 
one zero mode at each end of the chain. Thus the termi- 
nated end behaves like the ±1/0 interface and supports 
±e/2 excitations. A finite chain in the phase u = does 



not exhibit mid-gap states, indicating that the u = 
phase can indeed be identified with the vacuum. The 
corresponding DW states for model dU can be obtained 
from Eq. ([5]) by substituting 2t ^ A and 4Aa -^ vf- The 
— 1/0 DW behaves similarly to the 1/0 DW and does not 
require a separate treatment. 

In contrast to ±1/0 DWs, a —1/1 DW supports two 
zero modes. The numerical diagonalization reveals four 
mid-gap states (two per DW) approaching zero energy 
as the ring size is increased. The wavefunctions of these 
states are localized at the DWs. Analytical approxima- 
tion for the wavefunctions are obtained by linearizing (|2]) 
around k = 7r/2 and solving DW '0i/_i(x) and anti-DW 
'0_i/i(x) states 

V^1/_1(X) =5± e^^O ^^'^^^')a((2t5^K4A)-) |c±), (6) 

where |6±) = (| t)|^+) ^ | Dl^-)), \c±) = 
{\t)\y-)±\i)\y^)) . For -1/1 DW (1/ - 1) m{x) 
is increasing (decreasing) function satisfying m{xo) = 
and approaching constant value away from xq. Since 
there are two zero modes per DW (anti-DW), the 
possible charge states are, as in polyacetylene, ±e and 
depending whether both of the states are empty (e), one 
of them is excited (0) or both are excited (— e). In the 
case of an antisymmetric DM m{x — xq) = —m{xo — x) 
the valence band is spin paired so ±e charge states are 
actually spinless as in polyacetylene [8]. The argument 
for DW charge states is supported by calculations based 
on a method introduced in Ref. [12] recently. After 
an appropriate regularization, the DW charge (modulo 
multiples of e) follows from the expression 

q = —^ / dx duj dkQdkQ~^Qd^Q~^QdxQ~^ 

where Q = {iuj — H(k,x))~^ is the Wigner transformed 
Green function [9|, ll2|. To extract DW charge, we will 
use the continuum approximation of H in the vicinity of 
the DW. Evaluation of the DW charge formula conffims 
that possible charges states for a ±1/0 DW (—1/1) are 
half-integral (integral) multiples of e. 

The possible charge states of 1/0 and 1/ — 1 DWs coin- 
cide with spinless and spinfull polyacetylene DWs. How- 
ever, the corresponding wavefunctions are quite different 
from those of polyacetylene zero modes, where DW and 
anti-DW wavefunctions are nonzero only on separate sub- 
lattices [8] . In our model the sublattice and spin structure 
follows directly from the property C\ip^) = ±|'0^), where 
the different signs correspond to DW and anti-DW zero 
modes. Our numerical calculations indicate that mid-gap 
states are insensitive to even-odd variation in the num- 
ber of lattice sites and robust against weak disorder in 
system parameters. 

Physical realizations- Zigzag fermion chains are known 
to exist in various physical contexts. Structurally the 
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Figure 2: (a) Numerical solution of zero mode 
wavefunctions (blue triangles) for a ring with L = 1200 
sites, X/t = 1/8 and an analytical approximation based 

on continuum solutions (|5j). The red dashed line 

represents the magnetization profile, (b) In QSH wires 

a 1/0 DW supports a fractional charge state ±e/2 

distributed evenly between ±e/4 charges. 



most closely related examples are monoatomic Au, Pt, 
and Ir chains |13., J^]. For example, Au spontaneously 
forms zigzag chains with one s-type conduction electron 
per site. Since spin-orbit effects are pronounced in heavy 
elements, these metal chains could provide a realization 
of the studied phenomenon. The exact shape of these 
structures, magnitude of the spin-orbit hopping A and 
the ^-factor determining m depend on the details of the 
system and are not readily available. Thus it is difficult 
to give a realistic estimate for the energy gap setting 
the relevant temperature scale. However, if our model 
captures the essential physics for these systems, a finite 
chain exhibits ±e/2 edge states for any finite A and \m\ < 
2t at low temperatures. 

More concrete estimates can be given for the QSH 
nanowire model (j4]) which has similar topological proper- 
ties. In a thin wire the edge states on opposite edges have 
a weak overlap, inducing a gap A in the edge spectrum 
pi) . This gap depends sensitively on the width of the wire 
[16] and can be engineered by modulating it. Magnetic 
field can be employed to tune the magnetization m [15 1 
so by varying m and A it is possible to drive the system 
to different phases illustrated in Fig. [Tdl In Ref. [16] it 
was estimated that for wires of width 200 nm the inter- 
edge coupling is A = 5.22 K, while in Ref. [15] it was 
concluded that in-plane magnetization of m = 3.5 K can 
realized by magnetic fields of the order IT. These val- 
ues suggest that experimental investigation of the phase 
diagram and the zero mode excitations of QSH wires is 
possible within existing methods at sub-Kelvin tempera- 
tures and high magnetic fields. As discussed in Ref. [15|, 
the single-flavor Dirac spectrum realized in the edges of 
QSH systems have ±e/2 excitations located at magnetic 
DWs. In a QSH wire the ±1/0 phase boundary realizes 
another exotic possibility, bound pairs of ±e/4 charges on 



opposite sides of the wire, as illustrated in Fig. [20 The 
zero modes (J6j) are localized equally in both sublattices 
indicating that equal charges reside on opposite edges of 
the wire. Since the total charge on a DW is ±e/2, oppo- 
site edges must have equal charges ±e/4. The ±1/0 DW 
is taking place when magnetization and tunneling are 
equal, determining the critical width of the wire. In the 
narrow (wide) part of the wire tunneling can be made 
stronger (weaker) than the magnetization, resulting in 
the 0/1 DW when the wire width exceeds the critical 
width. Unlike e/2 charges studied in [15], e/4 charges 
rely on inter-edge coupling and are inherently bound to- 
gether. The maximum distance of e/4 charges is limited 
by the critical width of the wire where m exceeds A. 
Still, for wires of width 200 nm or wider the charges have 
a very small overlap and can be thought of as separate 
entities. 

Conclusion- We introduced a realistic ID fermion 
chain exhibiting three topological phases protected by 
chiral symmetry. Two of the phase boundaries support 
±e/2 zero modes while one of the phase boundaries have 
±e and neutral excitations. We showed that similar topo- 
logical properties are shared by quantum spin Hall wires 
in magnetic fields. We solved the zero energy states us- 
ing combination of analytical and numerical methods and 
discussed how predicted properties could be realized in 
atomic chains and quantum spin Hall wires. 
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